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Abstract: Neutralino dark matter in supersymmetric models is revisited in the presence 
of flavor violation in the soft supersymmetry breaking sector. We focus on flavor violation 
in the sleptonic sector and study the implications for the co-annihilation regions. Flavor 
violation is introduced by a single JIr — tr insertion in the slepton mass matrix. Limits on 
this insertion from BR(r — >> \i + 7) are weak in some regions of the parameter space where 
cancellations happen within the amplitudes. We look for overlaps in parameter space where 
both the co-annihilation condition as well as the cancellations within the amplitudes occur. 
In mSUGRA, such overlap regions are not existent, whereas they are present in models 
with non- universal Higgs boundary conditions (NUHM). The effect of flavor violation is 
two fold: (a) it shifts the co-annihilation regions towards lighter neutralino masses (b) 
the co-annihilation cross sections would be modified with the inclusion of flavor violating 
diagrams which can contribute significantly. Even if flavor violation is within the presently 
allowed limits, this is sufficient to modify the thermally averaged cross-sections by about 
(10-15)% in mSUGRA and (20-30)% in NUHM, depending on the parameter space. In 
the overlap regions, the flavor violating cross sections become comparable and in some 
cases even dominant to the flavor conserving ones. A comparative study of the channels is 
presented for mSUGRA and NUHM cases. 
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1 Introduction 

Supersymmetric standard models have a natural dark matter candidate namely, the lightest 
supersymmetric particle (LSP) if R-parity is conserved [1]. In mSUGRA/CMSSM models, 
the LSP typically is the lightest neutralino [2-4]. In most of mSUGRA /CMSSM parameter 
space, the lightest neutralino is mostly a bino (5°); the bino component being close to 
99%. With the bino cross-section being small, the neutralinos are overproduced resulting 
in a larger dark matter relic density compared to WMAP [5] allowed range. There are 
however, some special regions in the mSUGRA parameter space where the neutralino is 
able 1 to satisfy the relic density limits [7, 8]. These are the (i) Bulk region, (ii) Stop (t) co- 
annihilation region, (hi) Stau (f ) co-annihilation region, (iv) A— pole funnel region and (v) 

x See also Ref. [6] 
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Focus point/ Hyperbolic branch regions. The various processes which play an important 
role in each of these sub-cases is shown in Fig. (1). 
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Figure 1. Annihilation channels appearing in the 11m calculation. V and Z are the 

chargino and neutralino mixing matrices [9]. 

The stau-co-annihilation region requires the mass of the lightest stau, f\ to be close 
to the mass of the LSP. The stop-co-annihilation is typically realized with large A— terms, 
which is also the case with the bulk region [10]. Among the above depicted regions, dis- 
counting the case of large A— terms, f-co-annihilation and the focus point regions are 
most sensitive to pre-GUT scale effects and the see-saw mechanism [9, 11-15]. It has been 
shown that the co-annihilation region gets completely modified in the SU(5) GUT the- 
ory and leads to upper bounds in the neutralino masses [9]. Similarly, in the presence of 
type I, type II or type III see-saw mechanisms [9, 11, 14, 15] f-co-annihilation regions get 
completely modified. Strong implications can also be felt in the focus point regions unless 
the right handed neutrino masses are larger than the GUT scale [11]. GUT scale effects 
can even revive no-scale models [16]. It has also been shown that in the presence of large 
A— terms 'new' regions with f-co-annihilation appear [9, 17]. 

In the present work, we consider flavor violation in the sleptonic sector and study its 
implications for the co-annihilation regions. In generic MSSM, flavor violation can appear 
either in the left handed slepton sector (LL), right handed slepton sector (RR) or left-right 
mixing sector (LR/RL) of the sleptonic mass matrix. However, we concentrate on the 
flavor violation in RR sector as it has some interesting properties related to cancellations 
in the lepton flavor violating amplitudes as discussed below. Such flavor mixing is not 
difficult to imagine. It appears generically in most supersymmetric grand unified theories. 
A classic example is the SUSY SU(5) GUT model. If the supersymmetry breaking soft 
terms are considered universal at scales much above the gauge coupling unification scale 
(Mgut), typically the Planck scale, then the running of the soft terms between the Planck 
scale and the GUT scale could generate the RR flavor violating entries in the sleptonic 
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sector [18, 19]. 

For demonstration purposes, lets consider the superpotential of the SU(5) SUSY-GUT: 
W = hfjlOilOjSH + hfjlOiEjZH + • • • (i.i) 

where 10 contains {q,u c ,e c } and 5 contains {d c ,/}. As supersymmetry is broken above 
the GUT scale, the soft terms receive RG (renormalisation group) corrections between the 
high scale Mx and Mgut, which can be estimated using the leading log solution of the 
relevant RG equation. For example, the soft mass of 10 would receive corrections: 

- K)ro, - -T^2 A? Vu V tj (3m 2 , + Al) log U**-\ , (1.2) 

\ GUT / 

where Vij stands for the ij th element of the CKM matrix. Since 10 contains e c , the flavor 
violation in the CKM matrix (in the basis where charged leptons and down quarks are 
diagonal) now appears in the right handed slepton sector. Below the GUT scale, the 
RG scaling of the soft masses just follows the standard mSUGRA evolution and no further 
flavor violation is generated in the sleptonic sector in the absence of right handed neutrinos 
or any other seesaw mechanism. Assuming Mx ~ 10 18 GeV, the leading log estimates of 
the ratios of flavor violating entries to the flavor conserving ones, S RR = A^/m~, are 2 
given in the Table 1. We have taken Aq = and ht ~ 1. At 1-loop level 8 it is roughly 
independent of mo- 

Table 1. Flavor Violation generated in SU(5) Model 



1*1 


Value 




xRR 

Tjl 




7.8 • 10" 5 
2.0 • 10" 3 
1.4- 10 -2 



From the Table 1, we see that the RG generated 5^ R is typically of (9(10 -3 - 10" 5 ). 
Such small values will not have any implications on the co-annihilation regions or rare 
flavor violating decays. While non-universality at the GUT scale in this case is RG induced, 
there are models where non-universal soft terms can arise from non-trivial Kahler metrics 
in supergravity, this could be the case in models with flavor symmetry at the high scale a la 
Froggatt-Nielsen models (see for example, discussions in [20-26]). In such cases, the 5rrs 
could be much larger, even close to 0(1). These terms would then receive little corrections 
through RG as they are evolved from the GUT scale to the electroweak scale. Recently, 
in an interesting paper [27], supersymmetric models with Left-Right symmetry have been 
studied with particular emphasis on leptonic flavor violation. In these models, both left 
handed and right handed sleptonic sectors have flavor violation with the constraint that 

2 m~ is the flavor conserving average slepton mass. 
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^RR(^r) — ^LL(A r ), where A r is the left right symmetry breaking scale. In such cases it 
could be possible 3 to generate Srr ~ 0(1O _1 ). 

In this present work, we will follow a model-independent approach and assume the 
presence of a single flavor violating parameter A^? and study the implications of it for the 
co-annihilation region. We will consider the simplistic case of universal soft-masses at the 
Mqut scale with non-zero 5^ which is treated as a free parameter. To distinguish from the 
standard mSUGRA model, we will call this model 5-mSUGRA and similar nomenclature 
also holds for the other supersymmetry breaking models which we consider in this work. 

While flavor violating entries in the sleptonic mass matrices are strongly constrained in 
general, the constraints on leptonic 5^ entries are weak in some regions of the parameter 
space [29-31]. This leads to the possibility that large flavor violation could be present in 
the sleptonic right handed sector. In these regions cancellations happen between various 
contributions to the lepton flavor violating (LFV) amplitudes. If such cancellation regions 
overlap with regions where sleptonic co- annihilations are important, flavor violation has to 
be considered in evaluating the co-annihilation cross-sections in the early universe. This 
is the basic point of the paper where we show that flavor violating processes can play a 
dominant role in the co-annihilation regions of the supersymmetric breaking soft parameter 
space. The processes contributing to relic density in these regions are called flavored co- 
annihilations. 

It turns out that with mSUGRA/CMSSM boundary conditions, the parameter space 
where the flavor violating constraints are relaxed does not overlap with the f\ co-annihilation 
regions unless one considers extremely large values of 5 ^ 0.8. The overlap is not very sig- 
nificant and is mostly ruled out by other phenomenological constraints. However, if one 
relaxes the complete universality in the Higgs sector i.e., within non-universal Higgs mass 
models (NUHM), there is an overlap between these regions, paving way for large flavor 
violation to coexist with co-annihilation regions. 

The fact that in 5-NUHM these regions do overlap has already been observed indepen- 
dently by Hisano et al. [32, 33]. However, they have studied ji —> e^y transitions and their 
co-annihilating partner is not really a mixed flavor state. Further, they have not studied 
the relic density regions in detail. 

In this present work we elaborate on these regions and study the consequences of 
it. The rest of the paper is organized as follows: In section [2] we discuss the effect on 
5 in the co-annihilation regions both in the mass of the co-annihilating partner and in 
the cross section. We also show that overlap between regions of LFV cancellations and 
co- annihilations are not possible in 5-mSUGRA. In section [3] we show that in 5-NUHM 
regions do exist where flavored co- annihilations become important. Relative importance 
of various cross-sections in the flavored co-annihilation regions is elaborated in section 
[4]. We close with a summary and brief implications for LHC in [5]. In Appendix [A] we 
have written down the approximate expression of the soft-masses for mSUGRA and NUHM 
scenario for three different values to tan /3. In Appendix [B] we present 5-mSUGRA in more 

3 Subsequent to the appearance to this work on arXiv flavored co-annihilations have been studied by the 
group [28]. 
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detail using approximate results. Description of numerical packages used and numerical 
procedures followed are in Appendix [C]. In Appendix [D], we present loop functions which 
are relevant to the discussion in the text. In Appendix [E] we present the analytic form of 
the cross-sections for some scattering processes relevant for the present discussions. 




Figure 2. The Co-annihilation region with and without flavor mixing. In the above figure 
we plot the condition — m^o = for 5 = (blue line) and for 5 = 0.5 (green line). Here we 
have chosen tan /3 = 5 and Aq = 0. 



2 Co-annihilation with Flavor Violation 

Co-annihilations play an important role in reducing the (relic) number density of the dark 
matter particle by increasing its interactions at the decoupling point. It requires having 
another particle which is almost degenerate in mass with the dark matter particle and 
should share a quantum number with it [34]. In mSUGRA, Xi can have co-annihilations 
with f\ in regions of the parameter space where « m x°- We wn ^ now generalize this 
condition 4 in the presence of flavor violation. As discussed in the introduction, we will 
consider a single \i — r flavor mixing term in the RR sector, A^ to be present at the weak 
scale. Similar analysis also holds for the e — r flavor mixing. The slepton mass matrix is 
defined by 

dntD-l^Mf® (2.1) 
4 The condition can be more accurately expressed as = ra^o + 5m, where 5m lies within 10-15 GeV. 
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where $ T = |e L , (11, Tl, e R , p, R , tr| and 
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where, = rrif (Af — /itan/3)'s are the flavor conserving left-right mixing term, m 2 ~ 's 

are the left handed slepton mass term and m 2 ~ 's denote the right handed slepton masses. 

JR 

In the limit of vanishing electron mass 5 and zero flavor mixing in the selectron sector, 
we can consider the following reduced 4x4 mass matrix. This matrix is sufficient and 
convenient to understand most of the discussion in the paper. It is given by 



( 2 





m l L R 


\ 





m% 





tlr 


VLR 





VR 


Kr 


V 




RR 





(2.3) 



where, we have taken it to be real for simplicity. The lightest eigenvalue of the above matrix 
can be easily estimated. The lower 2x2 block can be diagonalized assuming that the flavor 
violating is much smaller than the flavor diagonal entries. A second diagonalization 
for the stau LR mixing entry can be done in a similar manner. This leads to a rough 
estimate of the lightest eigenvalue as: 



where 5 



^RR 



m~ R (l - S) - m r fi tan /3, (2.4) 
Requiring that the lightest eigenvalue not to be tachyonic, we find 



an upper bound on S as follows: 



S < 1- 



m r /xtan/3 



rn 



(2.5) 



TR 



This condition becomes important in regions of the parameter space where a ^> m~ and 

•R 

in regions where tan/3 is very large such that the second term approaches unity. For 
co-annihilations, 5 lowers the lightest eigenvalue of the sleptonic mass matrix. Non-zero 
5 shifts the 'standard regions' in mSUGRA towards lower values of M 1 / 2 , for a fixed 
rriQ. In other words, since the sleptons become lighter, the co-annihilations happen with 



5 In all our numerical calculations, we have used the full 6x6 mass matrix without any approximations. 
This approximation is valid only in models with universal scalar masses, like mSUGRA, NUHM etc. 
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lighter neutralino masses. To illustrate this point let us consider mSUGRA like universal 
boundary conditions at the GUT scale. The one exception to the universality of the 
scalar mass terms particularly slepton mass terms at GUT scale is in terms of the flavor 
violating mass term (A^ R ). We will call this model as 5-mSUGRA. Given that the 
parameter does not run significantly under RG corrections 6 , we can use the MSSM RGE 
with mSUGRA boundary conditions to study the low energy phenomenology. In Appendix 
[A.l], we have presented approximate solutions for the RGE of soft masses and couplings in 
mSUGRA. Using approximate formulae, in Fig. (2) we have plotted, the f— co-annihilation 
condition, m~o — ~ 0, with and without flavor mixing. We have chosen S = 0.0, 0.5 
and tan/3 = 5. As expected from the Eq.(2.4), the presence of flavor violating 8 shifts the 
co- annihilation regions more towards the diagonal in the mo — Mi plane. In table 2, we 

2 

show the spectrum for two points with 5 = and 6 = 0.5 which demonstrate that for fixed 
mo, a lighter neutralino can be degenerate with m^ in the presence of 5. 

Table 2. Spectrum in co- annihilation region with and without 5. 



Parameters Mass (GeV) 



mo 


200.0 


200.0 


Mi 

2 


1031.0 


458.0 


tan f3 


20 


20 


5 


0.0 


0.5 


m x\ 


439.22 


188.69 




439.24 


188.70 



Eq.(2.4) is a rough estimate and not valid for large 8. A more accurate expression 
is presented in Appendix [ B]. As we will see, this will not change the conclusions of the 
present discussion much. We will revisit this point again in the next section. 

The presence of 5 also affects the relic density computations in the co-annihilation 
regions . The thermally averaged cross section on which relic density crucially depends can 
get significantly modified with 5, where flavor violating scatterings are also now allowed. 
The typical f co-annihilation processes in the absence of flavor violation are XiXi ~~ ^ 
rf,/x/2,ee, Xi^i ^7, 'HTi — > tt, fif* rf, Xi^l ~~ ^ ^ T -> TiT* ~~ ^ 77- I n the presence 
of jiR — tr flavor mixing, the new vertices related to flavor mixing would contribute to 
the processes with flavor violating final states. The corresponding Feynman diagrams are 
shown in Fig. (3), where ji/r would mean that the final state could either be a ji or a r. The 
relevant Boltzmann equations for the neutralino and the lightest slepton (Zi), continue to 
remain as in the unflavored co-annihilation case, though the masses and the cross-sections 
appearing in them change. 

6 This is true as long as we stick to MSSM like particle spectrum and interactions. Additional interactions 
and particles can modify the flavor structure. 
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Figure 3. Co- annihilation channels appearing in the ft dm calculation with \i — r flavor violation 
in the right handed sector. Notice that there are now new final states where either \i or a r could 
appear. 

We have computed all the possible co-annihilation channels including flavor violation 
by adding the flavor violating couplings in the MSSM model file of well known relic density 
calculator, MicrOMEGAs [35]. The flavor violating co-annihilations contribute significantly 
to the total cross section and their relative importance increases with increasing 6 as ex- 
pected. So far we have not addressed the question whether such large flavor violating 
entries in the sleptonic mass matrix are compatible with the existing flavor violating con- 
straints from rare decay processes like r-^/i + 7orr^ jiee etc. Constraints from such 
processes have been discussed in several works. The constraints on right handed (RR) 
flavor violating sector are different compared to those of left handed (LL) sector as they 
only have neutralino contributions and have no chargino contributions. Furthermore the 
two neutralino contributions 7 can have cancellations amongst each other in certain regions 
of the parameter space as elaborated in refs. [29-31]. Following [30], the branching ratio 
for r —> n + 7 can be written as in the generalized mass insertion approximation 

These are the pure B° and the mixed B° - H° dia grams, as depicted in Fig. (4). 
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BR(r h. M7) = 5-78 x 10" 5 MkM^l x |^ (/b r _ Ir) \ 

iMl 



(2.6) 



where 7# 5j r and are loop functions are given in Appendix [D]. 

This amplitude is resultant from the two diagrams shown in the mass-insertion ap- 
proximation in Fig. (4). The first one is a pure Bino (B°) contribution whereas the second 
one is a mixed Bino-Higgsino (B — — H®) contribution. There is a relative sign dif- 
ference between these two contributions and thus leads to cancellations in some regions 
of the parameter space. In 5-mSUGRA, these cancellations occur when rrif R ~ 6 Mi or 
equivalently /i 2 ~ m~ R [30] . In regions outside the cancellation region the limit on Srr is of 
O(10 _1 ) for tan/3 = 10 and for a slepton mass of around 400 GeV [36] using the present 
on BR(r — » fi + 7) < 4.4 x 10 -8 [37]. In the cancellation region however the bound on 5 is 
very weak and S could be 0(1). 
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Figure 4. B° and B° — H° contribution in RR-insertion. The photon can be attached with 
the charged internal lines. 

A large 8 ~ (9(1) 8 would increase the flavor violating cross sections in the early uni- 
verse. The current bounds already push the value of 8 ~ 10 _1 for reasonable values of 
slepton mass ~ 400 GeV and tan/3 ~ 10. We look for regions where the bound is signifi- 
cantly weakened due to cancellations. This would require that there should be significant 
amount of cancellations among the flavor violating amplitudes to escape the bound from 
r — » n + 7- I n Fig. [5], we have presented the numerical results for mSUGRA with each 
panel representing a different value of 5 (0.2, 0.4, 0.6 and 0.8). tan/3 is fixed to be 20 and 
sign(/i) is positive. The details of the numerical procedures we have followed are presented 
in Appendix [C]. In all these plots, we have shown contours of BR(r /jlj) and the co- 
annihilation regions. The other constraints shown on the plot include, the purple region 
which is excluded as the LSP is charged, (m^ < ^ x o); the translucent black shaded region 
is excluded by search for a light neutral higgs boson at LEP, nth < 114.5 GeV, the light 
green region where the chargino mass is excluded by Tevatron, m ± < 103.5 GeV. The 
co-annihilation region has been computed including the flavor violating diagrams in the 
thermally averaged cross-sections. The relic density is fixed by the recent 7-year data of 

8 By definition S cannot be larger than 1. Here 0(1) means close to 1. 
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m [GeV] m [GeV] 




Figure 5. mo — Mi plane in 5-mSUGRA: The different contour shows branching ratio, BR(r — >• 
/jl 7) for 5 = 0.2, 0.4, 0.6 and 0.8 (from top left clockwise) and for tan j3 = 20, A = and sign(/i) > 0. 
The blue line indicates WMAP bound satisfied region. The black shaded region is excluded by direct 
search in LEP for the Higgs boson. The violet dots represent the present limits form LHC [38]. 
The red dot-dashed line indicates 1 TeV contour for gluino and blue dotted line marks the 1 TeV 
contours for first generation squark mass. The regions where the contours of BR(r — >• ^7) reaches 
< 10 -10 are the places where cancellations happen. In this region 5j^ R becomes unbounded because 
of the cancellation between the B° and B° — H° diagrams in Fig. (4). 



WMAP which sets it to be [5], 

n CD Mh 2 = 0.1109 ± 0.0056 (2.7) 

In the blue shaded region the neutralino relic density (Qdm) is within the 3a limit of [5], 
i.e., we require it to be 

0.09 < n DM h 2 < 0.12. (2.8) 

From the first panel of the figure, for 5 = 0.2 we see that there is no overlap in the regions 
where cancellation in the amplitudes for r — >► ^ + 7 happens (around BR(r — >► fij) < 10 -10 
) and the co-annihilation region (blue region). With increasing 5, as can be seen from 
subsequent panels, the co-annihilation region moves towards the diagonal of the plane as 
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the slept on mass becomes lighter, and the cancellation region which requires nrif R ~ 6 Mi 
also moves towards the diagonal with increasing 5. However, within 5-mSUGRA these two 
regions do not coincide except partially at the top end of the spectrum close to the upper 
bound of the the co-annihilation region. 





m [GeV] m [GeV] 



Figure 6. Panels (from top-right in clockwise direction) depict mo- Mi plane with raio = rrihd = 
0.5 • mo and m2o = mhu — 1-5 • mo for tan /? = 20, Ao = and sign(/i) > 0, with S = 0.2, 0.4, 0.6, 0.7 
respectively. Dark green region indicates inefficient REWSB. Purple region indicates l\ LSP. Black 
shade marks the region excluded by unsuccessful search by LEP, rrih < 114.5 GeV. The violet dots 
represent the present limits form LHC [38]. The red dot-dashed line indicates 1 TeV contour for 
gluino and blue dotted line marks the 1 TeV contours for first two generation squark mass. Blue 
strip bordering l\ LSP is the co-annihilation region. The different contour marks BR(r — >• jij). The 
regions where the contours of BR(r — >• /ij) < 10 -10 and below are the places where cancellations 
happen, which can be identified by their 'band' like structure. 

From Fig. (5) we can see that a very large (Sj^ R > 0.8) is required to make the cancel- 
lation region consistent with the co-annihilation region. In 5-mSUGRA having such large 
5 is consistent only very specific points of the parameter space (please see Appendix B for 
more discussion). Hence, we can infer from the above figures that within the 5-mSUGRA 
scenario the cancellation and co-annihilation region are disparate and no simultaneous so- 
lution exists. While the present discussion was based on numerical solutions for a particular 
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tan/3, one can easily convince oneself that it would be true for any tan/3 by looking at the 
analytical formulae. In fact, in the co-annihilation region, the branching fraction can be 
evaluated in the limit (m,f R — > Mi) and is given as 



BR(t -»■ fi-y) as 1.134 x 10 -6 x 



M^|cffif tan 2 /3 



(2.9) 



Mf 



where, we have used \a\ 2 « 0.5m? + 20M? and m% « m% + 2.5M-?. It is important 
to note that, the above expression obviously does not permit any cancellations. Thus 
within 5-mSUGRA, flavor violation in the co- annihilation region even if present would be 
constrained by the existing leptonic flavor violating constraints. In the following we see 
that this situation is no longer true in case, when, one relaxes the strict universality of the 
5-mSUGRA and considers simple extensions like non- universal Higgs mass models. 

Before proceeding to 5-NUHM, a couple of observations are important. Firstly, apart 
from the cancellation regions, the present limits on BR{r — >> \i + 7) constraint \b\ < 
0.11-0.12 for tan/3 of 20 and slepton mass of around 200 GeV (Mi ~ 500 GeV) in the co- 
annihilation regions. Since such values of 5 are allowed by the data, one can consider them 
to be present in 5-mSUGRA. A larger value of 5 would be valid for larger slepton masses. 
As discussed, this would lead to shifts in the parameter space of the co-annihilation region 
corresponding to mSUGRA. As a result, there is a shift in the spectrum also compared to 
mSUGRA. The thermally averaged cross-section are also modified. The shifts would be 
largest in the absence of any constraint from lepton flavor violation. For this reason, we 
look for overlapping regions between the cancellation and co- annihilation regions. Secondly, 
the cancellation region lies within a small narrow band. To the left and right of this band 
there could be regions of partial cancellations. These are present in Figs. (5). A discussion 
connected with this issue is present in Appendix [B]. 

3 Flavored Co-annihilation in 5-NUHM 

As we have seen in the previous section, in 5-mSUGRA, the fi parameter gets tied up with 
the neutralino mass in the co-annihilation region, thus leaving little room for cancellations 
within the flavor violating amplitudes. In the NUHM models, which are characterized by 
non- universal soft masses for the Higgs alone [39], the \i remains no longer restricted. This 
can be demonstrated with approximate formulae presented in Appendix [A. 2]. We denote 
the high scale mass parameters as m^(MouT) = m !o an( ^ m # d (^GUT) = m io- F° r tan/3 
= 20, using the approximate expressions in the Appendix [A. 2], we see that |/x| 2 has the 
form: 



Setting ml « m~ - 0.15M? and Mi w 0.41 lMi and taking the limit mf R -> Mi in the 

R 2 2 

co-annihilation region, we have 



fi\ 2 w 0.67 ml + 2.87 M? - 0.027 m 2 10 - 0.64 m| 



'20 



(3.1) 



H\ 2 ph 17 Ml - 0.027 m\ - 0.64 m| 



-20 



(3.2) 
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thus providing enough freedom 9 in terms of raio and ra2o to allow cancellations in the LFV 
amplitudes to co-exist with co-annihilation regions. 

The dark matter phenomenology of NUHM models has been studied by several authors 
[39-44]. The LSP is a neutralino in large regions of the parameter space and further, it 
can admit large Higgsino fractions in its composition unlike in mSUGRA. For simplicity, 
we concentrate on Bino dominated regions in the following. In such a case the lightest 
neutralino mass, in terms of SUSY parameters is as in mSUGRA: 

ra v o « 0.411Mi (3.3) 

Al 2 

For the lightest slepton mass one can use Eq.(2.4) where now m~ R at weak scale will 
be determined by the NUHM boundary conditions at the GUT scale. Similar to the 
mSUGRA case, approximate solutions can be derived for the NUHM case also and they 
are presented in Appendix (A. 2). Using the co-annihilation condition ra^ ~ ra % o and 
the cancellation condition m^ R « /x 2 , one can derive expressions for m 2 and ra^o where 
flavored co-annihilations are of maximal importance. 

The derived expressions for ra^ra^o are however, complicated. We found simpler 
parameterizations for regions where the LSP is Bino dominated and co-annihilations with 
the l\ are important. Examples of such regions are (i) 777,20 = 1.5 -mo and raio — 0.5 -rao and 
(ii) ra2o = 3 • rao and raio = ^o- For these values of raio and ra2o, flavored co- annihilations 
can exist for non-zero 5. In Fig. (6), we present in rao, Mi plane regions consistent with 
all constraints for 5 = 0.2,0.4,0.6 and 0.7, in an analogous fashion as to those presented 
in 5-mSUGRA section, Fig. (5). We have chosen ra2o = 1-5 • rao and raio — 0.5 • rao 
for this plots. The purple region is excluded as the LSP is charged, here mj < ra % o. 
Dark green region indicates no radiative electroweak symmetry breaking, |/x| 2 < 0. The 
translucent black shaded region is excluded by search for light neutral higgs boson at LEP, 
rrih < 114.5 GeV. As in 5-mSUGRA, we see that with increase in $23^5 ^1 — LSP region 
increases owing to the reduction of mass of l\. The impact of non- universality in the Higgs 
sector is negligible for ra^ in these regions. Analogously, regions excluded by light higgs 
search (ra^ < 114.5 GeV) are weakly affected in the presence of 5. Moreover, region with 
\fi\ 2 < is not affected by 5 as it is entirely governed by rao, raio and ra2o with maximum 
contribution from ra2o and rao- However, as expected the magnitude of BR(r — » ^7) 
governed by eq.(2.9), increases with The last panel of the figure shows regions where 
cancellation regions overlaps with the co-annihilation regions for 5 = 0.7. For a different 
set of values of raio and ra2o, for example, ra2o = 3 • rao, raio = mo the overlap regions 
can be found for even smaller values of 6. In these regions flavored co- annihilations play a 
dominant role. 

4 Channels 

The individual scattering processes involved in the computation of thermally averaged 
cross-section are called channels. The typical channels which are dominant in the co- 
annihilation region are hh — >> ZZ", Xih ~~ ^ ih Xih ~~ ^ Zh X1X1 ~~ ^ ^ W* ~~ ^ M etc - (they 
V| 2 « 20.5M! 2 in this limit in mSUGRA as can be seen from the expression below Eq. (2.9) 
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are about thirty of them in total). In the presence of flavor violation the number of 
these processes would be enlarged to include flavor violating final states. In the present 
section, we analyze the relative importance of the new flavor violating channels with the 
corresponding flavor conserving ones as a function of 5. To, do this we fix Mi and vary 

2 

5 and mo- In effect, this corresponds to the combination of horizontal sections of the co- 
annihilation regions of all the panels in Fig .(5) (Fig. (6)) for 5-mSUGRA (5-NUHM). In 
Fig. (7) we plot the dominant channels as a function of 8 in 5-mSUGRA. All the points 
satisfy relic density within WMAP 3a bound and lie in the co-annihilation region. Rest 
of the phenomenological constraints are also imposed, mo is varied from 100 to 600 GeV, 
whereas M 1 / 2 is fixed at 500 GeV, tan/3 = 20 and sign(/x) > 0. The Y-axis is percentage 
contribution to the thermally averaged cross section, (crv) defined by 

% (av)^ mn = ( 7 } ^ m " x 100 (4.1) 

(crv) total 



Table 3. —. — c anne for dominant channels for (5-mSUGRA 

{(TV) tota i 



Parameters 


Point I 


Point II 


Point III 




Mi = 500.0 GeV, 


Mi = 500.0 GeV, 


Mi = 500.0 GeV, 




tan P = 20, 


tan P = 20, 


tan P = 20, 




m = 165.6 GeV 


m = 169.6 GeV 


m = 249.0 GeV 


8 


0.197 


0.202 


0.5 




0.0910 


0.119 


0.120 




0.206 


0.227 


0.181 


Xi h -> 7 M 


6.53 x 10" 2 


7.47 x 10" 2 


0.13 


l\ l\ — y t t 


0.211 


0.181 


0.116 


l\ l\ — ¥ T fl 


0.130 


0.117 


0.165 


l\ l\ — > /i /i 


2.10 x 10" 2 


1.97 x 10" 2 


5.97 x 10" 2 


h ^77 


0.110 


9.65 x 10" 2 


9.93 x 10" 2 


X? h -»■ Z T 


5.67 x 10" 2 


6.23 x 10" 2 


4.96 x 10" 2 


Xi h^Zfi 


1.76 x 10" 2 


2.02 x 10" 2 


3.53 x 10" 2 


h l*^Zj 


5.00 x 10" 2 


4.42 x 10" 2 


5.18 x 10" 2 


Xi Xi -> t f 


2.02 x 10" 2 


2.81 x 10" 2 


2.27 x 10" 2 


XiXi^rft 


6.76 x 10" 3 


9.50 x 10" 3 


8.29 x 10" 3 


Xi Xi -> 1* P 


1.73 x 10" 2 


2.42 x 10" 2 


1.80 x 10" 2 



It should be noted that flavor violating constraints are not imposed for 5-mSUGRA in 
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• • • XxXx 


— >■ TT 


• * • XiXi 


-> Tjl 


• • • Xi Xi 





• • • 


h Xi - 


-> 7 T 


• • • 


h Xi - 


^ 7M 




Figure 7. Channels in 5-mSUGRA: The colored dots show relative contribution of a particular 
channel to ((Jv) to t' ^1/2 = 500 GeV and mo and 5 are varied to fit the co- annihilation condition. 
Here all the points satisfy WMAP 3<r bound (2.8). For the above plots tan/3 is fixed to 20 and 
sign(/i) > 0. Flavor violating constraints are not imposed here. 

this analysis. The current limits on BR(r —)*// + 7) constraint \6\ < 0.11 in the parameter 
space presented in the figure. For those values of 6 we see that the flavor violating channels 
contribute up to 5% of the dominant channel contribution. Larger values of 5 are not 



- 15 - 



allowed after the imposition of this constraint as there is no overlap between cancellation 
regions and co-annihilation regions in 5-mSUGRA. However, to study the features of the 
channels with respect to 5 it would be useful not to impose the BR(r — )> fi + 7) constraint 
for the present. 



Table 4. aV c ^ annel for dominant channels for (5-NUHM 

{(TV)total 



Parameters 


Point IV 


Point V 


Point VI 




Mi = 750.0 GeV, 


Mi = 750.0 GeV, 


Mi = 750.0 GeV, 




2 

tan p — zu, 


2 

tan p — zu, 


2 

tan p — zu, 




m = 199.3 GeV 


m = 216.0 GeV 


m = 592.1 GeV 


8 


0.01 


0.12 


0.767 




0.115 


0.116 


0.111 


Xi h 7 r 


0.190 


0.168 


0.116 


Xi h — >■ 7 A* 


4.74 x 10" 4 


3.89 x 10" 2 


9.89 x 10" 2 


Zi Zi — >• r r 


0.388 


0.280 


0.134 


h h t n 


1.90 x 10" 3 


0.127 


0.227 


h h H H 


2.39 x 10" 6 


1.48 x 10" 2 


9.37 x 10" 2 


h 


0.115 


0.123 


0.129 


X? fi ->• Z r 


5.50 x 10" 2 


4.88 x 10" 2 


3.35 x 10" 2 


Xi h Z [i 


2.02 x 10" 6 


1.11 x 10" 2 


2.28 x 10" 2 


h l$^Zj 


5.67 x 10" 2 


6.36 x 10" 2 


7.49 x 10" 2 


XiXi^rf 


1.14 x 10" 2 


1.13 x 10" 2 


3.72 x 10" 3 


Xi Xi -> t A 


2.80 x 10" 5 


1.77 x 10" 3 


3.53 x 10" 3 


Xi Xi -»■ A* A 


9.53 x IO- 3 


9.87 x 10" 3 


4.49 x 10" 3 



The upper left panel shows the % (av) for X1X1 ~^ ^? which contributes about < 5% 
total to (av) in this region of parameter space. In this case, the initial state masses are 
independent of 5 and mo, and thus, the only variation comes from the mass of intermediate 
state particle (l\). In Table (3), we presented the sample points which are represented in 
the plot. From the points, I and II of table 3, we see that a slight shift of 5 GeV in mo 
is still allowed by WMAP 3a limits, which changes the X1X1 cross-section by about 40%. 
This is the reason why the band of allowed points is broad in this channel. Other dominant 
channels are represented in subsequent panels of the figure. From the panel it is obvious 
that the dominant contribution comes from Xi h ~~ ^ 7 T an d h h t t channels. Each of 
which contribute to about 35% and 25% respectively to (av). Most of the flavor violating 
counterparts of these channels behave as expected, i.e. at large (5, they become comparable 
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Figure 8. Channels in 5-NUHM: The colored dots show relative contribution of a particular 
channel to (<Jv)tot- ^1/2 = 750 GeV and mo and 5 are varied to fit the co- annihilation condition. 
Here all the points satisfy WMAP 3<r bound (2.8). For the above plots tan/3 is fixed to 20 and 
sign(/i) > 0. Flavor violating constraints are imposed here, which causes the discontinuous regions 
in each of the channels. 

to the flavor conserving ones. One exception of this is the channel. Here the initial 
state composition crucially depends on '5' and also on tltr mixing. In such a situation, 



-17- 



its clear that the initial state cannot be attributed any flavor quantum number. In fact 
we find that the jl (smuon) component of l\ can be large ~ 50% even for 5 « 0.2 in some 
regions of parameter space. We see from the figure that the flavor violating final states 
dominates over the flavor conserving ones, as 5 grows beyond 5 ^ 0.2. The exact point 
of crossing of the flavor violating channels over flavor conserving ones is dependent on 
the parameter space chosen, crucially on tan/3 and fi. This is because the effective ftL^R 
and/or Jirtl coupling generated play an important role in determining the initial state 
composition. The last two panels shows some of the channels, which contribute negligibly 
to the (crv). In Appendix [E] we have given approximate formulae in m r /m /i limit for 
the dominant cross-sections. Using these and approximate formulae presented in appendix 
A features of full numerical analysis can be verified. More detailed analysis of cross-sections 
in the presence of flavor violation is various dark matter allowed regions will be presented 
elsewhere [45]. 



h Xi -> 7^ 
h Xi -> 7 M 
\ l\ — y t t 




b 



as 



10 1 - 



io u - 



io- 1 ^ 



0.74 



_i i i_ 



_i i i_ 



_i i i_ 



_i i i_ 



0.76 



0.78 



0.80 



0.82 



0.84 



c 23 



Figure 9. Dominant Channels contribution to the {crv) to t- Here tan/3 is fixed to 20, M 1 / 2 = 750 
GeV and mo and 5 are varied to fit the co-annihilation condition. Here all the points satisfy WMAP 
3cr bound (2.8). 

In Fig. 8, we present similar plots form channels in 5-NUHM case for the parametriza- 
tion chosen in the previous section. Here we have imposed BR(r — >> /ry) < 4.4 x 10 -8 to 
be satisfied along with relic density constraints. These channel show a similar pattern here 
as in 5-mSUGRA. However, as we can see from the panels, there is a gap between 6 = 0.2 



- 18 - 



to 5 = 0.7 where the parameter space does not satisfy BR(r — > fij) < 4.4 x 10 -8 . For 
points below 5 ^ 0.2, this constraint is satisfied as '5' is too small to generate appreciable 
r —> (i + 7 amplitudes. For 5 ^ 0.7, the constraint is now satisfied because of the overlap 
between the cancellation regions and co-annihilation regions. The relative contribution in 
the overlap region is magnified in Fig. 9 where all the channels contributions are pre- 
sented between 0.70 ^ 5 ^ 0.85. As we can see, flavor violating channels strongly compete 
with flavor conserving ones. A sample of the points in 5-NUHM is presented in Table 4, 
where points IV and V represent low 6 values whereas point VI represent the large 6 value 
signifying overlapping regions. 

Finally a note about relative contribution to relic density. We have 



tth 2 oc 



1 



[(Tv)tot J2 C 



av)i 

all channels 

l 



oc r -^- 4.2 

all canals 

For small 8 (~ O (lO -2 )) where (av) contribution to flavor violating channel is small, the 
estimate of relic density does not modify much from the flavor conserving case. However 
for large enough 5 (~ O (l0 -1 )), one tends to overestimate relic density, if one does not 
consider flavor violating scatterings while computing the thermally averaged cross-section. 



5 Summary and Outlook 

We have generalized the co-annihilation process by including flavor violation in the sleptonic 
\± — t (RR) sector. The amount of flavor violation admissible is constrained to be small 
by the limit on the BR(r — > fi + 7). This constraint is significantly weakened in regions 
of the parameter space where cancellations in the amplitudes takes place. We look for 
regions of parameter space where there is a significant overlap between cancellation regions 
and co- annihilation regions. The search is done in mSUGRA and NUHM augmented with 
one single flavor violating parameter in the \i — r (RR) sector. We found that while 
no significant overlap is possible in 5-mSUGRA, 5-NUHM allows for large regions where 
significant overlap is possible. 

The presence of flavor violation shifts the lightest slepton co-annihilation regions to- 
wards lighter neutralino masses compared to mSUGRA. While computing the thermally 
averaged cross-sections in the overlap regions, we found that flavor violating processes could 
contribute with equal strength and in some cases even dominantly compared to the flavor 
conserving ones. This is true even for 6 > 0.2 in some regions of the parameter space. 
Neglecting the flavor violating channels would lead to underestimating the cross section 
and thus in overestimating the relic density. A point to note is that if flavor violation is 
present even within the presently allowed limits, it could still change the dominant chan- 
nels by about 5% in 5-mSUGRA and more in 5-NUHM. Finally, We have probed only a 
minor region of the parameter space in the present work demonstrating the existence of 
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such regions. A comprehensive analysis of such regions and the associated phenomenology 
of their spectrum would be interesting in their own right. 

In this respect, a few comments on flavor violation at the LHC and ILC are in order. 
Detection of lepton flavor violation at the colliders like LHC is strongly constrained by 
experimental limits on rare lepton flavor violating decays. One standard technique to 
detect flavor violation at colliders is to study the slepton mass differences using end-point 
kinematics of cascade decays [46]. The typical sensitivity being discussed in the literature 

is i£?(h,lj) = ^== - O(0-l)% for e L - fa and 0(1)% for ji L - f L [47]. In the 

I y/ l i Ij 

presence of splittings are generated in all the three eigenvalues [48], e — , /i — r,e — r 
sectors. In the case discussed in this work, the typical splittings are (9(20)% to (9(70)% 
as the constraints from LFV experiments are evaded. Thus, far less sensitivity is required 
to measure these splittings compared to the regular case. Further investigations in this 
direction are however needed. Another interesting aspect of this scenario would be to 
measure widths for LFV decay processes like ~~ ^ Xi^flJ- These widths have been 
studied for the case of right handed slepton flavor violation in [49]. In NUHM, with a 
comparatively smaller value of \i one could expect large production cross sections for \\ 
and xf m the decays of colored particles. In fact, a full Monte Carlo study has been 
reported by Hisano et al. [33] for a particular parameter space point in the model. 

At the linear collider, it should be possible to identify the f co-annihilation region 
[50-54] by studying the polarization of the decay f\ — )> Xi T - I n the presence of flavored 
co-annihilations one should be able to see flavor violating decays of t\. Heavier particles 
like T2 and charginos would also have flavor violating decays. 

Finally lets note that we have considered the cancellations in the dipole operator of 
the t — y ji transitions, it does not guarantee us suppression in amplitudes associated 
with other operators. For example, in this region r— >/i7/ or r—^/ir/ could be sizable 
(~ 10 -9 — 10 -10 ) [55], which could be probed in future B-factories. Whereas, r — >> fij will 
continue to remain constrained and thus will not be detected. 

The focus of the present work has been to introduce new regions of parameter space 
where flavor effects in the co- annihilation regions could be important. More generally flavor 
effects could play a role in any dark matter 'regions' of the SUSY parameter space. Such 
studies are being explored in [45]. 
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A Approximate Solutions 
A.l mSUGRA Case 

In the approximation of small Yukawa couplings, we retain only Y tl I5, Y T and solve the 
RGEs semi- analytically. For the first two generations of the particles the dependence on 
tan f3 is very weak, so we take them to be valid for all tan /?. In deriving the approximate ex- 
pressions we have taken m t (Mz) — 165 GeV, m^Mz) = 3 GeV and m T [Mz) = 1-77 GeV. 
For tan/3 = 5, the first two generation masses at the weak scale are 

mg + 6.66M? (A.l) 

2 

m^ + 6.19M? (A.2) 

2 

ml + 6.22 Mi (A.3) 

2 

(m£)i )2 (M z ) ~ m^ + 0.51M? (A.4) 

2 

(m|)i )2 (M z ) ~ m^ + 0.17M? (A.5) 

2 

Third generation masses strongly depend on tan (3 than the first two generations. For low 
tan f3 = 5 their values are as follows 



(m£) 3 (M z ) 




- 0.036 Al + 0.65 + 0.16 A Mi + 5.66 M? 

2 2 


(A.6) 


(m 2 uh(M z ) 


rs_/ 


- 0.070 Aq + 0.31 ml + 0.30 A Mi + 4.26 M? 

2 2 


(A.7) 


(mlUM z ) 


rs_/ 


- 1.70 x 10 -3 Al + m| + 7.23 x 10 -3 A Mi +6.17M? 

2 2 


(A.8) 


(mlHMz) 


rs_/ 


- 7.34 x 10" 4 Al + ml + 6.29 x 10" 4 A Mi + 0.51 M? 

2 2 


(A.9) 


(mlHMz) 


rs_/ 


- 1.47 x 10 -3 Aq + mg + 1.26 x 10 -3 A Mi + 0.16 M? 

2 2 


(A.10) 




r-o 


- 3.30 x 10 -3 Aq + 0.99 ml + 0.01 A Mi + 0.48 M? 

2 2 


(A.ll) 


™ 2 h u (M z ) 




- 0.105 Ag - 0.046 + 0.46 A Mi - 2.95 M? 

2 2 


(A.12) 


\»\ 2 (M Z ) 




- 4158.72 + 0.110 A^ + 0.084 - 0.47 A Mi + 3.09 M? 

2 2 


(A.13) 



(m^)i, 2 (M z ) ~ 
(m 2 D \ 2 {M z ) ~ 
(m^)i )2 (M z ) ~ 
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For medium tan f3 = 20 their values are as follows 



(m 2 Q ) 3 (M z ) ~ 


- 0.048 A 2 , + 0.62 ml + 0.20 A Mi + 5.54 M? 

2 2 


(A.14) 


(m^) 3 (M z ) ~ 


- 0.070 A 2 , + 0.33 ml + 0.30 4 Mi + 4.32 M? 

2 2 


(A.15) 


(m£) 3 (M z ) ~ 


- 0.023 A 2 , + 0.91 m 2 , + 0.10 4 Mi + 5.86 M? 

2 2 


(A.16) 


(m|) 3 (M z ) ~ 


- 0.011 Al + 0.97 ml + 8.38 x 10- 3 A Mi +0.50M? 

2 2 


(A.17) 


(m|) 3 (M z ) ~ 


- 0.021 Al + 0.93 m 2 , + 0.017 A Mi + 0.15 M? 

2 2 


(A.18) 


^L( M ^) - 


- 0.046 A 2 , + 0.83 ml + 0.16 A Mi + 0.01 M? 

2 2 


(A.19) 




- 0.105 A 2 , - 0.007 ml + 0.46 i Mi - 2.86 M? 

2 2 


(A.20) 


l/f (M z ) = 


- 4158.72 + 0.106 A 2 , + 0.009 m 2 , - 0.46 A Mi + 2.87 M? 

2 2 


(A.21) 



For high tan (3 — 35 their values are as follows 



(m 2 Q ) 3 (M z ) - 


0.058 A 2 , + 0.53 ml + 0.25 A Mi + 5.26 Af 2 

2 2 


(A.22) 


(r4) 3 (M z )~- 


0.064 Al + 0.33 m 2 , + 0.27 A Mi + 4.35 M ? 

2 2 


(A.23) 


(m 2 D ) 3 (M z )c- 


0.052A^ + 0.727m^ + 0.23A Mi +5.26M? 

2 2 


(A.24) 


(m|) 3 (M z ) ~ - 


0.027 A 2 , + 0.89 m 2 , + 0.02 A Mi + 0.49 Af ? 

2 2 


(A.25) 


(m|) 3 (M z )c- 


0.055 A 2 , + 0.78 m 2 , + 0.03 A Mi + 0.12 Af 2 

2 2 


(A.26) 




0.105 A 2 , + 0.48 m 2 , + 0.36 A Mi - 0.91 M? 

2 2 


(A.27) 


m 2 Hu (M z ) ~- 0.095 A 2 


- 0.005 m 2 , + 0.41 A Mi - 2.81 M? 

2 2 


(A.28) 


|/i| 2 (M z ) = - 4158.72 + 0.095 A 2 , + 0.005 m 2 , - 0.41 A Mi + 2.81 Af ? 

2 2 


(A.29) 



A.2 NUHM case 

In our notation m\o — mH d {MouT) and m2o = rnH u {MGUT)- For tan/3 = 5, at the weak 
scale the first two generation masses are 

(m 2 Q ) h2 (M z ) ~ Wq + 6.66 Af 2 + 0.009 (mf - m| ) (A.30) 

2 

(m 2 D ) 1:2 (M z ) ~ mo + 6.19 M? + 0.018 (m 2 — m| ) (A.31) 

2 

(m^)i j2 (M z ) ^ + 6.22 M? - 0.036 (m? - m^ ) (A.32) 

2 

(m|)i )2 (M z ) ~ m-o + 0.51 Af? - 0.027 (m 2 - m| ) (A.33) 

2 

(m|)i >2 (M z ) ~ m-o + 0.17 Af? + 0.053 (mf - m| ) (A.34) 

2 

Third generation masses strongly depend on tan f3 than the first two generations. For low 
tan f3 — 5 their values are as follows 

(m 2 Q ) s (M z ) ~ - 0.036 + 0.77mg + 0.16A Mi + 5.66 M? +7.90 x 10~ 3 m 2 10 

2 2 

- 0.125 mlo (A.35) 
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(mfj) 3 (M z ) ~ - 0.070^ + 0.54 m§ + 0.30 4)Mi + 4.26M? - 0.035 mj 

2 2 

- 0.196 m^o (A.36) 
(m 2 D ) 3 (M z ) ~ - 1.70 x 10 _3 A§ + mg + 7.23 x 10 _3 A Mi + 6.17 M? + 0.016 m? 

2 2 

- 0.018 m^ (A.37) 
(m|) 3 (M z ) ~ - 7.34 x 10" 4 ^+ m 2 , + 6.29 x 10" 4 AoMi +0.51M? - 0.027m? 

2 2 

+ 0.027 m 2 2Q (A.38) 
(m|) 3 (M z ) ~ - 1.47 x 10 _3 A§ + ml + 1.26 x 10" 3 AoMi + 0.16 M? + 0.052 m 2 

2 2 

- 0.053 m| (A.39) 
m ff d ( M z) - - 3.30 x 10 -3 Al - 7.32 x 10 -3 m 2 , + 0.01 A Mi + 0.48 Af ? + 0.969 mf 

2 2 

+ 0.027 m| (A.40) 
m^ (M z ) ~- 0.105 A 2 , -0.70 mg + 0.46 A Mi — 2.95 M\ + 0.027 mf 

u 2 2 

+ 0.625 m^o (A.41) 
|^| 2 (M Z ) = -4158.72 + 0.110 A 2 , + 0.72m 2 , -0.47 A Mi + 3.09 M? + 0.012 m 2 

2 2 

- 0.650 m^ (A.42) 
For medium tan /3 = 20 their values are as follows 

(m 2 ) 3 (M z ) ^ - 0.048 Ag + 0.75 m 2 , + 0.20 A Mi + 5.54 M? - 6.30 x 10 -3 m 2 w 

2 2 

-0.120ml (A.43) 
(m 2 r ) 3 (M z ) ~- 0.070 Al + 0.55 m 2 , + 0.30 A Mi + 4.32M? - 0.034 m? 

2 2 

-0.190m^ (A.44) 
(m 2 D ) 3 (M z ) ~- 0.023 Ag + 0.94 m 2 , + 0.10 A Mi + 5.86M? - 0.015 m 2 w 

2 2 

- 0.015 m^o (A.45) 



l io 



(m|) 3 (M z ) ~ - 0.011 A 2 , + 0.987^ + 8.38 x 10 _3 A Mi + 0.50M? - 0.038m? 

2 2 

+ 0.027 m| (A.46) 
(m|) 3 (M z ) ~ - 0.021 A^ + 0.95 m 2 , + 0.017 A Mi + 0.15 M\ + 0.030 m 2 

2 2 

- 0.053 mlo (A.47) 
m 2 H (M z ) ~ - 0.046 Ag - 0.11 mo + 0.16 A Mi + 0.01M? + 0.913 m 2 w 

2 2 

+ 0.030 mlo ( A - 48 ) 
m 2 Hu (M z ) ~-0.105Ag-0.67m^ + 0.46A Mi -2.86M? + 0.030 mf 

2 2 

+ 0.634 mlo ( A49 ) 
H 2 (M Z ) =- 4158.72 + 0.106 A 2 , + 0.67 m 2 , - 0.46 A Mi +2.87 M? - 0.027 mf 

2 2 

-0.636 mlo ( A5 °) 
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For high tan f3 = 35 their values are as follows 



{m 2 Q UM z ) 



0.058 A 2 + 0.69 ml + 0.25 A Mi + 5.26 Af ? - 0.037 mf 



0.120 m| 



(ml) 3 (M z ) ~- 0.064^ + 0.55m§ + 0.27 + 4.35M? - 0.029mf 



0.194 m| 



(m^) 3 (M z ) ~ - 0.052 Al + 0.82 m§ + 0.23 A Q Mi + 5.26M? - 0.081 mj 



(mlUMz) 
(mlUM z ) 



0.010 m^ 

0.027 Aq + 0.93 + 0.02 A Mi + 0.49 Af? - 0.063 m\ 



+ 0.027 rr4o 



0.055 + 0.85 mg + 0.03 A Mi +0.12 Aff - 0.019 mf 

2 2 

0.054 m|o 



m 2 H (M z ) a- 0.105 Ag-0.35 + 0.36 A Mi — 0.91 M? + 0.789 m 



10 



+ 0.038 mi 

m 2 Hu (M z ) - 0.095 A 2 - 0.67 + 0.41 A Mi - 2.81 Af? + 0.036 m 



(A.51) 
(A.52) 
(A.53) 
(A.54) 
(A.55) 

(A.56) 



10 



+ 0.629 ra^o 



(A.57) 



M 2 (Mz) 



4158.72 + 0.095 Ag + 0.67 ml - 0.41 A Mi + 2.81 Aff - 0.036 m| 



0.629 m| 



(A.58) 



B Lightest Slepton Mass in <5-mSUGRA at Large 6 



From the plots presented in section 2, Fig.( 5), for the case of (5-mSUGRA, the following 
two things can be inferred: (a) the co-annihilation condition increasingly moves towards 
the diagonal in (mo,M 1 y 2 ) plane with increasing 5 and (b) the cancellation region are 
almost independent of the value of 8 in (mo, M 1 / 2 ) plane. The question then arises if there 
is some region at large '5' where the two regions coincide. In the present appendix, we 
explore this question. The analysis presented here is based on the approximate solutions 
of Appendix [A.l] and we will comment on the full numerical solutions at the end of the 
section. 

The effective 4x4 matrix of eq. (2.2) can be diagonalized as follows. First the lower 
2x2 block is rotated by an angle 6, given by, 

2A RR 



tan 261 



A" 2 2 ' 

m~ — mi 



(B.l) 



The eigenvlaues of this lower block can be easily read off from the mass matrix. They are 



A± ~ 2 



2 + 4A^ R 



(B.2) 
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Figure 10. 5 Contours: Upper bounds on S in various of the parameter space using the non- 
tachyonic condition. 



For ~ m f R (which is true for low tan/3 regions), the eigenvalues have the following 



form: 



\ 2 ± ~rh 2 ± A RR 



rh\l±5 RR ) 



(B.3) 
(B.4) 



where m 2 = - (m? + m 2 R ). Next we have to diagonalize the tlr entry. The eigenvalues 



after this rotation are approximately given as 



(m% + A 2 ) ± A /(m? L - A 2 _)2 + 4 cos 2 9^ A? 



» T 'LR 



(B.5) 



In the limit (m 2 L — A 2 ,) ^> Af LR (the corresponding angle is very small in this limit) 10 , 
which is the case for large 5, we can write the above eigenvalues as 



o o o o f 2 COS 2 UT A ~ 

« + a 2 _) ± « - x_) i + K r A 2_;r 



(B.6) 



So, the lightest eigenvalue of the effective 4x4 mass matrix of eq.(2.2) is given as 

cos 2 IIT A? 



ri ~ a 



J » T <LR 



m% - A 2 



COS 2 Our A % 



m\ L - m 2 (l - S RR ) 



(B.7) 
(B.8) 



- 25 - 




200 400 600 800 1000 
m [GeV] 




1 -I 

200 400 600 800 1000 

m [GeV] 



Figure 11. Cancellation and Co-annihilation region in <5-mSUGRA 



Which essentially suppresses the left-right mixing term compared to eq. (2.4). And de- 
manding the lightest eigenvalue to be non-tachyonic we get an upper bound on 5rr as 
below 

cos 2 9 UT A? 

Srr<1 / 2 TLR (B.9) 

Which matches with eq. (2.5) in the limit cos 9^ —> 1. 

In Fig. (10) we have plotted the tachyonic condition (R.H.S of eq. (B.9)) using the 
approximate results of Appendix [A]. It has been plotted for two values of tan/3 20 and 
35. The contours represents the upper bounds on S in those regions of the parameter space 
to avoid tachyonic leptons. As we can see, increasing tan/3, tightens the bound a bit. 
In Fig. (11) we have shown the cancellation condition /i 2 ~ m f R an d the co- annihilation 
condition ~ m~o for two values of S = 0.8 and 0.9. In both the pannels, the brown and 
magenta solid lines indicate co- annihilation condition for 5 = 0.8 and 0.9 respectively. The 
green dashed line satisfy the cancellation condition, whereas the orange and red dashed 
lines satisfy the cancellation condition with the n parameter being 30% corrected than 
its tree level value. Comparing the Figs. (10) and (11) we can see that there could be 
some points which could evade both the tachyonic condition as well as have cancellations 
amongst the LFV amplitudes and still satisfy the co-annihilation condition. However in 
practice in full numerical calculation, we could not find any points consistent with both 
these conditions as other phenomenological constraints rule them out. As can be seen 
from the figure, a 30% correction to the /i parameter could shift the overlapping region to 
very small values of (mo,M 1 / 2 ) or no overlap at all for 8 ~ 0.9. This approximates the 
implications of adding the full 1-loop effective corrections to the SUSY scalar potential. 
However, the co-annihilation region could allow for partial calculations in LFV amplitudes. 
Such regions are difficult to distinguish in a numerical analysis. 
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We will now return to Eq. (B.5) and consider the limit (rrv~ — A 2 ,) <C Af M , which is 
an interesting limit as it is relevant for the regions which appear in channels plots discussed 
in section 4. From eq. (B.6), there could be a value of 5 as well as parameter space in 
(mo, tan/3) where ra~ ~ \ 2 _. In these regions, the corresponding mixing angle is 

very large and the subsequent diagonalization is very different. It turns out that at least 
three mixing angles in the slepton mass matrix are large in this parameter space. The plots 
presented in Figs. (7) and (8) contain these regions. More details of these regions will be 
discussed in [45]. 

C Numerical Procedures 
C.l SuSeFLAV and MicrOMEGAs 

The numerical analysis is done using publicly available package MicrOMEGAs [35] and 
SuSeFLAV [56]. SuSeFLAV is a fortran package which computes the supersymmetric spec- 
trum by considering lepton flavor violation. The program solves complete MSSM RGEs 
with complete 3x3 flavor mixing at 2-loop level and full one loop threshold corrections 
[57] to all MSSM parameters and relevant SM parameters, with conserved R-parity. Also, 
the program computes branching ratios and decay rates for rare flavor violating processes 
such as fi — )> e7, r — >► e7, r — » [i 7, fi — » e7, \i~ — >► e + e~ e~, r~ — > fi + \i~ fi~ , r~ —> e + 
e~ e~, B —± 57 and (g — 2) M . 

In the present analysis we use M? ole = 173.2 GeV, M v h ° le = 4.23 GeV and M? ole = 
1.77 GeV. In determining the lightest higgs mass (ra/J we use approximations for one loop 
correction which are mostly top-stop enhanced [58]. We use complete 6x6 slepton mass 
matrix to correctly evaluate the inter-generational mixings and masses in the presence of 
flavor violation. 

Moreover we consider flavor violating couplings stemming from lepton flavor violation 
in the RR sector of f — fl. 

(a) (b) 

k y% 



Ix Xa Ix Xb 

Figure 12. (a) Neutralino-slepton-lepton vertex and (b) Slepton- lepton-chargino vertex. 

• Neutralino-slepton-lepton: 

The interaction Lagrangian for neutralino-slepton-lepton is written as 

£ = k (Xf AX P L + Y,f AX P R ) xa Ix + h.c. (C.l) 
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Where the coefficients are defined as 

^iax = K i [cos9 w (0 N ) A 2 + sm0 w (O N ) A i] U X ^M W cos (3 - m h cos9 w (O n ) A3 U x ^ 3 

(C.2) 

and 

Zf AX = -K x [2sin<9 w M w cos/3 U x ,i+s (O n ) a1 + m k cos9 w U x ,i (O n ) A3 ] (C.3) 
where 

e 1 
Kl y/2 sin 9 W M w cos /3 cos 9 W 

The Interaction Lagrangian for chargino-slepton-neutrino is 

C = Pi (nf BX P L + Ylf BX P R ) x + b lx + h.c. (C.5) 
Where the coefficients are 

nfsx = -^V (0 L ) B iU x ,i (C.6) 

Sin C7w 

smf/w v2Mwcosp 

Where, J7x,z is the 6x6 matrix which diagonalizes the sleptonic mass matrix, here 
the indices i = 1 to 3 and X = 1 to 6. (Oat) Am is the 4x4 neutralino mixing 
matrix, where A, m = 1 to 4 and (Ol) Bn is the 2x2 chargino left eigenvector matrix, 
where B,n = 1, 2. mj. is the mass of the lepton l{. In our notation Pl — and 

These couplings are programmed into MicrOMEGAs through CalcHEP [59] package. 

C.2 Constraints Imposed 

• We check for efficient radiative electroweak symmetry breaking, requiring |/i| 2 > 
for valid points. 

• We require rrif > m x o as LSP is neutral. Regions for which this condition is not true 
is excluded as f LSP regions. 

• We impose lower bounds on various sparticle masses that results from collider exper- 
iments. m h > 114.1(GeV), m x ± > 103.5(GeV) and m f > 90(GeV) [60]. 

• 2.0 x 10- 4 < BR{b 57) < 4.5 x 10~ 4 [37]. 

• We also check for the D-flat directions, while checking for the EWSB condition and 
charge and color breaking minima [61-63]. 
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D Loop Functions 



In this appendix we define the relevant loop functions that contribute to the amplitudes of 
flavor violating leptonic process BR(r — » ^7) as presented in the appendix of [30] 



XL 



Ml 



rn 



2 ' 



xr 



Ml 



rn 



2 > 

R 



VL 



m 



2 ' 



VR 



A* 



rn 



R 



Ib r and Ir are defined as follows, 



lBM M h m h m R 



2 2 



VR hi(x R ) 



VlQ\{x l ) - y R gi(x R ) 



I R (m z R ,Mlfi 2 ) 

m* R y R - x R 
The functions g\ and h\ are defined as follows, 

1 - x 2 + 2xln(x 



1 



VR 



[hi{x R ) - hi(y R )} 



E Cross-Sections 



-, hi(x) = 



1 + Ax - 5x 2 + (2x 2 + Ax) ln(x) 



(D.l) 



(D.2) 
(D.3) 

(D.4) 



In this appendix we present the approximate formulae for the relevant cross sections. We 
do not attempt to discuss a complete comparison of the analytical expressions and full 
numerical results in the present paper, that is left for an upcoming publication. These 
expressions generalize the existing expressions [64] in the literature to include full flavor 
violation in the sleptonic sector. The expressions are presented only for the dominant 
channels and in the limit m r ,m /i — » 0. More detailed expressions and their simplifications 
will be discussed elsewhere [45]. 



Numerator 

O channel 7^ : 7 V-^'-U 

Denominator 



E.l WW -> TT 



The cross-section of l\ l\ —¥ r r process is as follows. This process involves t- and ^-channel 
Xi exchange. In the following and in rest of the cross-sections e is the electric charge, #w 
is the weak mixing angle and A% is the mass of the W-boson. We get the simplified form 
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of the above cross-section in the limit of m r — >> as below, where the numerator is 



S 2 + S 2 _Ws (s - 4m?J - Y^ti (s + 2m| ; - 2m 



x log 



s + 2m 2 - 2m? - \ sis — 4m? 



5 + 2m 2 - 2m? + * sis — 4m? 



5 + 2m? - 2m? 



log 



s + 2m 2 - 2m? — \ s Is — 4m? 



s + 2m| ; - 2m? + Js (.-4m?) 



x <! 2£^m~o + 2X^X 2 _m^ 



+ m| ; + S^) 2 5 - 4£^m? 



5 { 5 — 4m? 



2 ( min + m~ + m? ( s — 2m? 



4E 2 E 2 m|o 



+ 2£ 2 £ 2 _m 2 5 - 8£t£lm z ~o m~ 

' Xi ' Xi «i 



^2 vh2 2 



And the denominator is 

327rs M^f cos 4 f3 cos 4 #w sin 4 #w (s — ^ m f^j 
Following appendix C the coupling structure is: 
. f-x°i-h: 

K X (E+Pfl + E_P L ) 



where 



and 



S + = [cos#w ON {1,2) + sin0 w CW(1, 1)] cos/3 M w Z7(l,3) 



£_ = -2sin<9 w M w cos/3 17(1,6) CW(1, 1) 
Where K\ is already defined in eq. (C.4). 



(E.2) 



(E.3) 



(E.4) 



(E.5) 



(E.6) 
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E.2 hh fir 

The simplified form of the l\ l\ fi r cross-section in the limit of m r , m M is calculated 
below. This process involves t- and ^-channel Xi exchange. The numerator of the cross- 
section is 



(E 2 , A 2 + E 2 _A 2 ) yjs (s-4m?) - (E 2 A 2 + E 2 A 2 ) 



x (s + 2m|o - 2m? ) log 



5 + 2m 2 - 2m~ - A 


Is ( S -4m?) 


s + 2m? - 2m~ + A 
Xi h \ 


/ s (s — 4m~ ^ 



5 + 2m? - 2m~ 

Xi /i 



log 



s + 2m? 2m~ - \ s (s — 4ra~ 
Xi h V V «i 



S + 2m| ? - 2m? + ^ ( S -4m?) 
(S 2 , A 2 + S 2 A 2 , ) m 4 ? + (E 2 A 2 + E 2 A 2 ) m{ 
+ m|o ((E 2 , + S 2 _) (A 2 + A 2 , ) s - 2 (S 2 , A 2 _ + S 2 A 2 + ) m?) } 

Js(s- 4m 2 ) f 

^ _ { 2 (E 2 A 2 + E 2 A 2 ) m 4 

+ mf + m 2 is -2m 2 )) I Xl 



( mi n + mi + m? n ( s — 2m? 
V Xi h Xi V Zi, 

+ 2 (E 2 A 2 + E 2 A 2 ) m£ + (E 2 (2A 2 - A 2 ) - E 2 (A 2 - 2A 2 .)) m| ? s 



} 



- 4 (E 2 A 2 + E 2 A 2 ) m| ? m? 

And the denominator is 

32tts M w cos 4 /3 cos 4 6> w sin 4 6> w (s - 4m? ) 
Here the coupling structure is: 

K x (A + P R + A_P L ) 



(E.7) 



(E.8) 



where 



and 



A+ = [cos0 w CW(l,2) + sin0 w ON(l, 1)] cos/3 M w 17(1,2) 



(E.9) 



(E.10) 



2sin0 w M w cos /3 17(1,5) CW(1,1) 



(E.11) 
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E.3 XiXi ~+ r/flr 

In the limit m T —> the cross-section for Xi Xi ~~ ^ t t is calculated. This process involves 
t- and ^-channel l± exchange. The numerator is 



( S -4m| ; ) 



sm l + { m h- m l) 



-{ (£ 4 + 4E 2 + Si) sm? + 2 (£ 4 + 3E 2 £ 2 _ + Ei) 



-5 + 2m? - 2m~ 

Xi /i 



log 



"V 5 ( 5 - 4 ^lo) + (^-2m|o+2m?) 



s ( s — 4m|o ) + ( s — ^ m |5 2m [ 



|s (-2E 2 £ 2 m |o + (E 4 + 4E^ E 2 + E 4 ) m 2j + 2 ^4 + 3E 2 s ^ + s 4 ) 



And the denominator is 

128vr s cos 4 j5 sin 4 6\y cos 4 9 W (s - 4m| ? ) 



(E.13) 



Whereas in the limit rn T . m M — >► the cross-section for Xi Xi ~~ ^ M T i s calculated. This 
process involves £- and ^-channel l± exchange. The numerator is 



e 4 < 



(s - 4m| ? ) 



sm~ + 



T I {£ 2 (2A 2 _ + A 2 ,) + E 2 (A 2 + 2A 2 + )} sm\ 



+ {E 2 . (3A 2 _ + 2A 2 ) + E 2 (2A 2 _ + 3A 2 ,)} (m? ? - m?)' | 



-5 + 2m? 2m~ 

Xi /i 



log 



s (s - 4m|o) + (s - 2m?o + 2m? 



S ( S -4m| ? ) + ( S -2m 2 ; + 2m?) 
{E 2 (3A 2 _ + 2A 2 + ) + E 2 (2A 2 _ + 3A 2 ) } (m 2 ; - m?)' 

+ s {- (E 2 A 2 + E 2 A 2 ) m| ? + (E 2 (2A 2 _ + A 2 ) + E 2 (A 2 _ + 2A 2 + )) m; 

And the denominator is 

128tt s cos 4 (3 sin 4 W cos 4 % (s - 4m| ; ) 



(E.14) 



(E.15) 
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■4 Xih -»■ 7t//* 



E 

In the limit m T — >■ the cross-section for %5 ^1 — 7 r is calculated. This process involves 
s-channel r mediation and t-channel l\ exchange. The numerator is 



(^ + £ 2 _)e 4 < 



log 



>|? + ml) + y + (-8 + m?) 2 - 2m 2 ; (s + m?) _ 

4o - 3m? ) + (s - 2m? + 2m? ) Wml + (s + m? ) 2 - 2m 2 (s + m? 



And the denominator is 

327T 5 cos 2 /3 sin 2 #w c °s 2 #w | m ^o + ^— 



Whereas in the limit m M — » the cross-section for x? 'i ^ 7 M is calculated, 
involves s-channel ^ mediation and t-channel /i exchange. The numerator is 



(E.16) 

x? ( s + m l)} 
(E.17) 

r ° 7 - ^ " is calculated. This process 



s + m~ 



• 2m : 



(A^ + A 2 _)e 4 log 



m 



o- (s + ra? ) - \/m~o + 



5 + m~ 



S (m|o-3m?) + ( S -2m|o + 



2m? 



mlo + {—s + w-fj) 2 — 2m|o (s + m~ 



(E.18) 



And the denominator is 



327T s cos 2 f3 sin 2 #w cos 2 #w |^~o + (^—s + m~Q — 2m| (s + m~Q | 

(E.19) 
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